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Entropy barriers and ageing states appear in martensitic structural-transition models, slowly re¬ 
equilibrating after temperature quenches, under Monte Carlo dynamics. Concepts from protein 
folding and ageing harmonic oscillators turn out to be useful in understanding these nonequilibrium 
evolutions. We show how the athermal, non-activated delay time for seeded parent-phase austenite 
to convert to product-phase martensite, arises from an identihed entropy barrier in Fourier space. In 
an ageing state of low Monte Carlo acceptances, the strain structure factor makes const ant-energy 
searches for rare pathways, to enter a Brillouin zone ‘golf hole’ enclosing negative energy states, 
and to suddenly release entropically trapped stresses. In this context, a stress-dependent effective 
temperature can be dehned, that re-equilibrates to the quenched bath temperature. 

PACS numbers: 64.70.Q-, 81.30.Kf, 64.70.K-, 87.15.Cc 


I. INTRODUCTION 

Ageing states of glassy systems are a longstanding 
pnzzl^. They are pictured as being metastably trapped 
in a multi-valley free-energy landscape in configuration 
space. Their re-equilibration times {t ^ gAF/T| 
inverse rates) to find the global minimum, depend on 
the set of free energy barriers {AF = AU — TAS} be¬ 
tween basins. Although both energy barriers AU and 
entropy barriers |A5'| = —AS can contribute, these 
two delay sources are distinct. Energy barriers require 
thermal activation, with Arrhenius rates as Boltzmann 
factors t~^ ^ , while entropy barriers do not 

have such inverse temperatures in the exponent. Hence 
if deep quenches face unsurmount able energy barriers, 
t ^ >> then re-equilibrations could be dom¬ 
inated by the entropy barrier delays t ^ coming 

from searches at almost constant energy AU/T « 1, 
for rare canyon-like pathways connecting the basins. 

The folding of proteins involve such entropic searches. 
The folded state of the macromolecule is one of a huge 
number of configurations, and a random search would 
take a very long time. Like a golf hole in a fiat surface, 
there are many ways of failing, and only a few of succeed¬ 
ing: this is the entropy barrier. The observed rapid fold¬ 
ing is understood through the concept of a guiding funnel 
leading to the folded stat^^l®. Protein folding toy models 
of a Brownian particle searching outside a golf hole (‘un¬ 
folded state’) for a funnel inside it (‘folded state’), also 
find entropy barriers, at the golf hole edge^. 

Entropy barriers and ageing behaviour, arise in several 
non-interacting statistical models. Backgammon models 
of freely hopping particleJ^, and ageing harmonic oscilla¬ 
tors under Monte Carlo dynamic^, both without energy 
barriers, can show glassy behaviour from entropy bar¬ 
riers alone. In an ageing state of harmonic oscillators, 
the MC aeeeptanee fraetions in a sweep over all sites 


are very small, and decrease monotonically and slowly 
with time. Kinetically constrained models with trivial 
statics, but with imposed constraint on the spin-flip dy¬ 
namics, can exhibit slow relaxations, and glass-like freez¬ 
ing at ‘entropy catastrophes^^. Glass-like quasi-steady 
states are found in models with infinite-range or power 
law interaction^. 

We find that models of structural transitions, used to 
study unusual quench responses of athermal martensites, 
can also shed light on such general issues of nonequilib¬ 
rium statistical mechanics. 

Martensitic transitions can have unusual dynamics; 
have naturally non-uniform states and power law in¬ 
teractions; and can be a rich source of interesting 
nonequilibrium modelJ^H^^. These structural transitions 
have strain-tensor components as the order parame¬ 
ters, with parent-phase high-symmetry ‘austenite’ con¬ 
verting to different variants of lower-symmetry product- 
phase ‘martensite’, separated by elastic domain walls. 
These twin-boundaries are oriented in preferred crystallo¬ 
graphic directions by anisotropic power law interactions, 
arising from generic St Venant Compatibility constraints, 
that ensure lattice integrit}^!!. 

Martensitic materials are classified as ‘isothermal’, 
with activated, slow conversions; or ‘athermal’, with very 
rapid conversions when quenched below a martensite 
start temperature, and with no conversions, above iP. 
However in puzzling experiments, some athermal materi¬ 
als convert even above the rapid-conversion temperature, 
after delays of thousands of second^. 

Monte Carlo (MC) simulations in athermal parame¬ 
ter regimes, of a square-to-rectangle (or a 2D version 
of the tetragonal-to-orthorhombic) transition, have con¬ 
sidered re-equilibration after a temperature quenclP^, 
with continuous strains represente d as i n various con¬ 
texts by discrete-strain pseudo-spin d^^ l ^^ i For a square¬ 
cell austenite converting to one of two possible rectangle- 
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cell variants of martensite, the pseudospin can take on 
three values of 0, dzl, as in a Blume-Capel modeJ^. 

For athermal-regime parameterJI^, dilute martensitic 
seeds in austenite of initial fraction n^(0) << 1, are sys¬ 
tematically quenched to well below the scaled Landau 
temperature T = Tq = 1. They quickly form a marten¬ 
sitic droplet, that searches for an autocatalytic twinning 
channel^, and finds it after an austenite-martensite con¬ 
version delay , when there is a sharp rise of the marten¬ 
site fraction towards unity. The conversion delays can 
be very long or very short, depending on teniperature. 
The simulations^ correctly model experimentS^, show¬ 
ing very fast (slow) conversions at low (high) T. In this 
paper, we provide an understanding of these martensitic 
re-equilibrations, that are also very relevant for glassy 
systems. 

The conversion times are insensitive to energy-barrier 
scales, and therefore can only arisS^^ from entropy bar¬ 
riers tjn ^ What is the nature of the entropy 

barrier that blocks immediate access to states of man¬ 
ifestly lower energy? How can the barrier crossings at 
nearby temperatures, be very fast, or very slow? 

We answer these questions, using concepts from pro¬ 
tein folding and ageing harmonic oscillators. The natu¬ 
ral description is in Fourier space. The initially isotropic 
strain strueture faetor makes a delay-inducing search on 
a constant-energy surface, for an anisotropic ‘golf hole’ 
in the Brillouin zone. The golf hole is bounded by the 
/c-space line where the energy spectrum of martensitic 
textures vanishes. We identify the distortion pathways 
for the structure factor to cross the entropy barrier and 
enter the negative-energy funnel region, when there is a 
sudden spike in the MC acceptances, with trapped stress 
released as heat. A strain-related effective temperature 
can be defined. 

In Section II we outline the pseudospin model used 
previously. In Section III we analyse the textural evolu¬ 
tion in Fourier space, and in Section IV understand the 
crossing of the vapour-liquid entropy barrier. In Section 
V we consider domain-wall thermodynamics, with details 
in the Appendix. Finally Section VI has a summary and 
conclusions. 


II. THE PSEUDOSPIN MODEL 

The strain pseudospiiP^^^ Hamiltonian is derived from 
a free energy F{e) in the order parameter strains e(r), 
evaluated at the Landau minima through e e{T)S{f). 
Here s{T) is the strain magnitude at the minima; and 
the pseudospin locates the triple-well minima as S{f) = 
0,±1. Thus F{eS) = H{S) = Hl F Hq F He, where 
the Landau term is Hl ^ ^gL{T)S‘^{r); the Ginzburg 
term is Hq ^ and the St Venant Compat¬ 

ibility term is He ^ {Ai/2)^V{r — r)S{r)S{r). The 
interaction is an anisotropic power law in the separation 
R that is scale-free, V{XR) ^ V{R)/X^^ with a spatial 


average that is zero. Here A is a difference operator on 
the square reference lattice. 

In Fourier space A ^ iK where = 2sin(/c^/2). 
The Hamiltonian is diagonal in Fourier space, and is for¬ 
mally that of ^-labelled inhomogeneous oscillator d^ l ^^ l 
The Hamiltonian energy of a pseudospin texture from 
MC dynamics at a time t is 

= (2.1a) 

k 

with a dimensionless martensitic-strain spectrum 

e{k) ^ 9l{T) + ilK^ + T(1 - <5,-o)C(fc), (2.16) 

where the (positive) interaction kernel i^^^ 

V{k) = {Kl - Klfl[K^ F {^A,/A^){K,Kyf]. (2.1c) 

The kernel V ^ cos^ 20 = (1 + cos4^)/2 vanishes along 
favoured 0 = ±7r/4 or = Fky diagonal directions, and 
so a nonzero contribution He r\j T,yik)\S{k,t)\^ 7 ^ 0 is 
a domain-wall mz5-orientation energy. The Fourier kernel 
average over the Brillouin zone (BZ) is V = [V] 0.3. 

Her^^Ai,A 3 are elastic constants; Eq is an elastic 
energy per unit cell; and Kq{T) = 2Eoe‘^{T)^ where 
s{T) = [|{1 + The Hamiltonian en¬ 

ergy for the k = 0 uniform state is just the Landau 
energy H = NEofLY.^S‘^{r), where /l = and 

griT) = r — l + (c^ — 1)^ < 0, favouring martensite 
for r < 1. Here the scaled temperature 

r(T) = (T-Te)/(To-Te) (2.2) 

at the first-order Landau transition temperature Tq is 
unity r(To) = 1; while at the Landau spinodal Tc where 
metastablel^ austenite becomes unstable r(Tc) = 0. No¬ 
tice that H depends on the quenched bath temperature, 
as a bulk pseudo-spin can flip to equilibrate to T, in a 
single t ^ t F 1 increment. The Hamiltonian H{t) then 
describes the slower pattern evolutions of domain walls 
separating the S = FI variants; or separating one variant 
and the S = 0 austenite. In the glass terminology, the do¬ 
main wall patterns are ‘inherent structures’ at landscape 
minimgP^. Similar models coupling strain to magnetisa¬ 
tion or to random disorder would be relevant for other 
glassy phenomensP. 

Powerlaw isotropie interactions ^ 1/R^ with a < d, 
with divergent spatial averages, show unusual statistical 
behaviour, with non-extensive thermodynamic functions; 
inequivalence of ensembles; and ergodicity breaking of 
quasi-stationary states whose lifetimes can diverge with 
system siz^. The a = d ‘marginal’ case has logarith¬ 
mically divergent averages. In our case, the Compatibil¬ 
ity interaction between the order parameters arises from 
minimizing non-order parameter strains sub ject to a con¬ 
straint connecting derivatives of all strain^^^^^. Thus 
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the ^ = 0 interaction contribution is zero as in ( 2 . 1 c). 
The spatial average then vanishes, ^ 

(1 — %o)%o ~ ^ 5 t^ 6 -niarginal case. 

Nonetheless, we find finite incubation delays from finite 
entropy barriers. 



T 


FIG. 1. Temperature-Time-Transformation: Log-linear plots 
of mean martensitic conversion times tm{T) and domain-wall 
orientation times tc{T) versus quench temperatures T. Pic¬ 
tures are snapshots of domain-wall ’vapour’ {tm > t), ’liq¬ 
uid’ (tc > t > tm) and ’crystal’ (t > tc)- Quench regions 
are indicated: Region 4 with T 4 > T > T 2 ; Region 2 with 
T 2 > T > Ti; and Region 1 with Ti > T > Tg. 

Re-equilibration under a quench-and-hold protocoP^, 
is described by the dynamic structure factor 

p{k,t) = \S{k,t)\‘^, (2.3) 

and its BZ average is the martensite fraction: 

^p{k,t)/N = J2sHr,t)/N = n^it). (2.4) 

k r 

Fig 1 shows schematic curves from the datsP^^^, of 
the Temperature-Time-Transformation (TTT) evolution 
of the seeded system, quenched to a temperature T.The 
free-running system is monitored as it passes through 
domain wall phases of vapour, liquid and crystal. The 
martensitic conversion time where = 0.5, is 

also the TTT phase boundary between vapour and liquid. 
There is also another time tc{T) for the orientation of 
domain walls, that will be considered elsewhere. 

We will consider three temperatures in three distinct 
temperature RegionP^. For T 4 > T = 0.76 > T 2 in Re¬ 
gion 4, the re-equilibration is dominated by the vapour- 
to-liquid or conversion dela}®! t^(T). For T 2 > T = 
0.55 > Ti in Region 2, the total delay also has contri¬ 
butions from the liquid-to-crystal delay tc{T). Finally, 
for Ti > T = 0.4 in Region 1, the conversion time is 
negligible, and the domain-wall orientation time tc{T) 


dominates. In Region 3, even though the quenched tem¬ 
perature is still below the Landau transition temperature 
To = 1 > T > T 4 = 0.824, the small, dilute martensitic 
seeds disappear into surrounding metastable austenite, 
and are not re-nucleated. Hence in Region 4 we study a 
T = 0.76 sufficiently below T 4 , for a reasonable number 
of the Nruns = 100 lo convert, in a reasonable time. 


III. TEXTURAL EVOLUTION IN FOURIER 
SPACE 

Fig 2a) shows that for a quench to T = 0.76 the Hamil¬ 
tonian energy of ( 2 . 1 ) is nearly flat, with H{t) ^ 0, up 
to t ^ tsm- It starts to go negative at t ^ ti ^ 0.85 tm] 
and falls rapidly at t The fall slows down at 

t ^ tsci and finally flattens at t ^ tc- By contrast, 
for a quench to T = 0.4 in Region 1, the energy drops 
almost immediately. Fig 2b) shows the austenite frac¬ 
tion no(t) = 1 — nm{t) behaves similarly, incubating at 
no{t) = no(0) 1 for T = 0.76 before falling; while being 

expelled immediately, for T = 0.4. 



FIG. 2. Energy and austenite fraetion evolutions: For T = 
0.76,0.4 in quenched-temperature Regions 4 and 1 respec¬ 
tively, the falloff of a) the Hamiltonian energy, H{t)/N vs t; 
and b) the austenite fraction no(t) = 1 — nm{t) vs t. 


Since the evolving energy is first zero, and then neg¬ 
ative, the relevant spectrum is thus a zero-energy plane 
through the e{k) surface, plus the negative energies be¬ 
low it. Fig 3 shows that the resulting relief plot nat¬ 
urally depicts a flat surface containing a golf hole de¬ 
fined by e{k) = 0, and a funnel e{k) < 0 inside it. The 
sidebar shows the temperature-dependent, anisotropic 
golf hole edge, that is large (small) at low (high) tem¬ 
peratures. There is an outer (inner) squared-radius of 
- 2/Re{T) ± 6 with 6 = {AiV/2^^) > 0 , 
where the average is 2/Rc = —(fli/Co) — <5 > 0 . 

Protein folding is understood through concepts such as 
golf holes and funnels in configuration spaced. Here, we 
find such concepts appearing naturally in martensitic re¬ 
equilibration, but in a more easily represented form in the 
label space k of strain modes, in which the Hamiltonian 
is diagonal. For protein folding, this would correspond to 
the label space of folding normal-modes, of the protein 
model Hamiltonian. 
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FIG. 3. (Color online) Golf holes and funnels. Relief plot for 
T = 0.76 of relevant martensitic spectrum e{k) versus k in the 
Brillouin zone, with a zero-energy plane. Sidebar: Golf hole 
edges e{k) = 0, for quench temperatures top to bottom, that 
are low T = 0.42 (explosive conversions); medium T = 0.76 < 
T 4 (slow, rarer conversions); and high T = 0.91>T4 = 0.824 
(no conversions, divergent entropy barriers). 



FIG. 4. (Color online) Strueture faetors for domain-wall 
phases: Relief plots of ln[l + \S{k,t)\‘^] versus {kx,ky) in the 
Brilliouin Zone for a) vapour; b) liquid; c) crystal. Insets: 
Corresponding contour plots. 


Fig 4 shows the Fourier-space structure factors 
{p[k,t) = \S{k,t)\‘^} as relief plots, for domain-wall 
phases of ‘vapour’, ‘liquid’, and ‘crystal’. Fig 5 shows 
the MC evolutions that transform the phases into each 
other. 

We first discuss the insets of Fig 5, that are the evolv¬ 
ing coordinate space textures {5'(r, t)} as previousl}ff^, 
but now labelled by the characteristic times of Fig 2. 
See the movie in Supplementary Material. The in¬ 
sets show that the random seeds quickly form a vapour 
droplet of zero energy, that fluctuates in place up to a 
time t = tsm ^ 100- it is in an ‘incubating’ or age¬ 
ing stat d^ l ^ l ^^ . The single-variant droplet then finds 
and enters an energy-lowering, autocatalytic-twinning 
channeP^of alternating opposite variants, around a time 
t = ti 220 . At a time t = tm ^ 240, a domain-wall 
liquid forms, with walls of wandering orientation. After 
a symmetry-breaking choice of a diagonal at t = tsc ^ 
260, a domain wall crystal of oriented twins forms, be¬ 
yond t = tc ^ 350. 

The main Fig 5 shows the same evolving textures, but 
now in Fourier space. See also Fig 4 and Fig 6. For 


t the ageing state p{k^ t) for the incubating vapour 

droplet persists unchanged as a broad, isotropic gaussian, 
poised over the butterfly-shaped golf hole of Fig 3. One 
might expect the gaussian to promptly distort along di¬ 
agonals, to fit into the correspondingly anisotropic golf 
hole, and narrow, to enter the negative energy funnel. 
Surprisingly, it does not do this. For t > tsm r\j 100 , 
it waits in an incubation stage, to develop wings along 
kx^ ky axes. See also the contour plots of Fig 6, at these 
times. Then for t < ti the peak narrows and then rises 
sharply, and for tm > t > G, enters the golf hole, where 
it adopts the bi-diagonal symmetry of the funnel. After 
a symmetry-breaking at t ^ tsc^ the structure factor at 
t > tc IS along a single diagonal. Thus the Fourier space 
distribution develops mis-oriented wings along the axes, 
before it forms wings along the Compatibility-favoured 
diagonals. In coordinate space, the austenitic S = 0 
spins at the droplet surface must flip collectively to pro¬ 
duce S' = ±1 surface spin regions of the right symmetry: 
much like a kinetic constraint!^, but here self-generated. 
The improbability of finding this collective-spin distor¬ 
tion constitutes the entropy barrier. 

Fig 6 shows the contour plots in the BZ corresponding 
to the relief plots of the main Fig 5. The p{k, t) value of 
a point on such contours represents the Fourier intensity 
or ‘occupancy’, at a given k. It would also be interesting 
to monitor the evolving occupancy at a given energy. We 
define the energy occupancy distribution /)(e, t), similarly 
to that of a protein-folding simulatioiP, 






Fig 7 shows the evolving histogram of the single-run 
p(e, t) versus energy e. The arrow marks the energy e = 0 
of the golf hole edge. The negative energies 0 > e > gp 
are the funnel region. The distribution remains fixed, 
up to t = tsm = 100, and then a small positive energy 
peak appears, corresponding to the wings along the axes 
of Fig 6. (It appears at different onset times in different 
runs, so a time average would wash it out.) The weight 
of the distribution moves more into the negative energy 
region, and by around t = tm = 240, it is almost entirely 
in the funnel. Note the long-lived, occupancy ^ike at 
the energy of the golf hole edge and its environ^. This 
disappears at on crossing of the entropy barrier. 

As shown in Fig 8 the final distribution, for all 
quenches, is an inverse-energy falloff in the excitation en¬ 
ergy above the bulk Landau term, e = e — gp >0: 


pie,t-,T) ^ l/e. (3.2) 


For a continuous-variable displacement of a harmonic os¬ 
cillator, averaging with a Boltzmann factor would yield 
the same inverse-energy behaviour. The form is also inde¬ 
pendent of Hamiltonian energy scales Eq and anisotropic 
stiffness constants Ai. Since the domain walls are sparse, 
and discrete, the energies are also discrete. The energy 
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FIG. 5. (Color online) Textural evolutions: Snapshots of Fourier-space strain structure factors as ln[l + t)|^] for successful 

conversions, after a quench to T = 0.76. Insets show coordinate-space textures {S{f,t)}, with green (blue/red) denoting 
austenite (martensite variants). See the movie. The martensite fraction conversion time is tm{T). Other crossover times 
tsm,ti,tsc,tc are defined in the text, and Fig 2. 



t=200 ti=220 t^240 t^=350 


FIG. 6. Contour plots of strueture faetor in Fourier spaee: 
Colour contour version of Fig 4 snapshots of ln[l + \S(k, t)|^] 
for T = 0.76. The three domain wall phases: vapour, liquid, 
crystal correspond to three contours: isotropic gaussian (e.g 
t = 80), X-shaped (e.g t = 200), single diagonal or fan-shaped 
(e.g t = tc — 350). 


has an upper cut off, that is estimated^ as log^^Q (e) < 
9.1, consistent with the simulation results. 

Notice that although the distribution p{k, t) or p(e, t) 
of Fig 5 or Fig 7 has much of its weight poised above 
the funnel states inside the golf hole, these domain-wall 
modes labelled by k or e, do not immediately collapse into 
available negative energy states. Re-equilibration does 
not follow a strategy of ’every mode for itself’. Rather 
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FIG. 7. Time evolution of single-run energy oeeupaney dis¬ 
tributions : Plots of p{e,t) versus e for quench to T = 0.76 
for a single run, at various times. Arrow marks the golf hole, 
of energy e = 0. The funnel region is 0 > e > = ~\9l\- 

A small peak at large energies appears on the formation of 
the entropically critical droplet for t > tsm — 100; and an 
occupancy spike from an entropy barrier persists at the golf 
hole energy right up to t = tm — 240. 


there is an ‘all modes together’ strategy: the modes first 
partially equilibrate so there is no net inter-mode en¬ 
ergy exchange, setting up some nonequilibrium mode- 
distribution; followed by a slower emergence, as entropy 
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FIG. 8. Final energy oeeupaney distributions: Log-log plot of 
p(e, t) versus e = c—ql showing final 1/e behaviour, regardless 
of quench temperature, and of different energy scales Eq = 
3,4, 5, 6. Solid line has a slope —1.0; and an estimated upper 
energy cutoflP^is ~ 9.1, consistent with simulations. 


barriers are crossed, of an equilibrium mode-distribution 
at the bath temperature. 

IV. CONVERSION DELAYS: ENTROPY 
BARRIER CROSSING 

We need to understand how conversion delays from en¬ 
tropy barriers, can be so drastically different, at nearby 
temperatures. At low temperatures of Region 1 of 
T < Ti = 0.42, seeds in austenite convert explosively 
to martensite variants, for every run. At very high tem¬ 
peratures of Region 3 of Tq = 1 > T > T 4 = 0.82 there 
is a complete ‘blocking’ of conversion to energy-lowering 
martensite, with the zero-energy seeds dissolving for ev¬ 
ery run, into zero-energy austenite. And for Region 4 of 
T 4 > T > T 2 = 0.68 there is a rise in both the conversion 
time and the fraction of blocked runs on approaching T 4 , 
so the mean time diverges like a Vogel-Fulcher lavff^ 

An understanding of entropy barriers |A5'| that are ei¬ 
ther zero, sharply rising, or infinite, comes from quenches 
to Region 4, with a Fourier distribution p{k, t) that starts 
as an isotropic gaussian and ends as an inverted fan along 
one diagonal, as in Figs 4, 5. We parametrize the distri¬ 
bution through a weight 77 of the emerging anisotropy. 

The distribution is separated, outside the golf hole (t < 
ti), and inside the funnel (t > ti), as 

p{k, t) = pg{k, -t)+ pf{k, 7])e{t - ti). (4.1) 

We will focus on the constrained search pathways outside 
the golf hole for t < ti, through the parametrization 

Pg = V{n^(0)(l+6^7/(t))+77(t) cos46>}C^e“^^/^'^^ (4.2) 


where Cg normalises the gaussian to unity, and bg = 
^m(0)~ . The evolution parameter 77 (t) then 

carries the fourfold anisotropy of the kernel of ( 2 . 1 c), as 
pg ^ 77 (t)(l+ cos4^). Here the normalisation (2.4) yields 
p{t) = [(n^(t)/n^(0)) — mhg > 0. The distribution is 
isotropic with p{t) = 0 , during the n^(t) = n^( 0 ) incu¬ 
bation for tsm > t- Whereas for ti > t > tsm ^ nonzero 
?7 > 0 induces an angular modulation, that increases the 
distribution at 0 = 0 , 7 r/ 2 .., i.e. along the ky axes. 

Writing the /c 7 ^ 0 Hamiltonian energy of (2.1) as av¬ 
erages < ... > over the distribution (4.2), so 

H/Nn^ = 42 <K^> +gL + (Ai/2) < V{k) >, (4.3a) 

we obtain on the zero-energy plane outside the golf hole, 
a constraint linking the Ginzburg, Landau and St Venant 
contributions, 

H/Nn^eo = {2<j'^-2/Rc}+{v{t)/nm{t)}5/2 = 0. (4.36) 

In the 77 = 0 ageing state the average golf hole radius 
determines the gaussian width or inverse droplet size as 
2 ( 7 ^ = 2/Rc- From the constraint of (4.3b), any decrease 
in width must be compensated by an increase in mis- 
orientation energy ^ 5 ^ AiV oi the last term: the 
77 > 0 wings must indeed, first emerge along the axes, 
before the diagonals. This explains the observed p(/c, t) 
distortions, of Figs 5,6 for ti > t > tsm- At t = G, 
when 77 (^ 1 ) = 77 ^(ti), the width from (4.3b) narrows to 
2 ( 7 ^ = {2 /Re) — b/2] and then to the inner radius. If 77 is 
(unphysically) taken to be negative, favouring diagonal 
wings right away, then the constraint of (4.3b) makes the 
width larger, going in the wrong direction. 

The delays tm ^ are understood through the 

temperature-dependent golf holes in the sidebar of Fig 
3. For T below Ti, the golf hole in the BZ is large, and 
the flat distribution from seeds directly forms a liquid 
distribution of Fig 4, entering the funnel immediately 
for every run. The conversion time is negligible, and 
its entropy barrier is zero, |A5'| 0. For T approach¬ 

ing T 4 the golf hole shrinks, and hence the search times 
rise; the entropy barrier diverges as |A5'| ^ 1/(^4 — T), 
and the fraction of runs converting to martensite falls, 
yielding Vogel-Fulcher behaviour. For T > T 4 the golf 
hole inner radius G?_ closes, and the resulting 4-petalled 
golf hole topology presents the isotropic gaussian with an 
infinite entropy barrier. Thus even though the marten¬ 
site Landau energy is lower than the austenite energy for 
To > T > T 4 , it becomes ergodically inaccessible to the 
small and dilute initial seeds. 

For nucleation by activation over energy barriers, a di¬ 
vergent droplet timescale is associated with a divergent 
size in coordinate space. By contrast, for non-activated 
entropy barrier crossing, a divergent search timescale for 
droplet pathways is associated with a shrinking bottleneck 
in Fourier space. 

Under MC dynamics for a given run 77 , the total Hamil¬ 
tonian energy H(t) of a texture {5'(r, t)} goes to H{t-\-l) 
at the next MC sweep. The probability for a given energy 
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change /S.E to occur at a time t, from Hamiltonian incre¬ 
ments AH{t,n) = H{t + 1) — 77(t), is obtained through 
an average over all runs n = 1, 2...Nruns' 

^ Nrur^s 

P{AE,t) = — - ^ ^AH(t,n),AE' (4.4a) 

runs 

n=l 

At early times, the probability is peaked at negative val¬ 
ues AE = — |A^| < 0, with an asymmetric shoulder 
on the negative sid^ and at long times, this becomes 
an equilibrium distribution, symmetric around zero (not 
shown). To determine if there is a dominant energy 
change during the evolution, regardless of when it oc¬ 
curs, we average the energy release probability over the 
entire holding time, P{AE) = ^^p\AE^t)/t^. 

Fig 9 shows P{AE) for various temperatures and en¬ 
ergy reductions. There are a few large magnitudes of 
energy release, but mostly, P{AE) falls as a powerlaw 
in the magnitudes ^ l/\AEp, with a common expo¬ 
nent 7 2. This suggests that the domain-wall adjust¬ 

ments have no characteristic energy scale, and are like 
small earthquakes, of all scales. Acoustic emissions oc¬ 
cur in martensites, from twin boundaries inducing energy 
changes, and power law distributions have been seerP^ 
with exponent close to 2.3. 
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FIG. 9. Probability distribution of energy releases: Log- 
log probability of energy releases at any time, P(AE), versus 
magnitude of releases |AF|. 


The Monte Carlo acceptance fractioiP Aacc(0: is the 
fraction of N sites where the MC move is accepted during 
a given sweep at t. Ageing non-interacting oscillators, 
have small and monotonically decreasing A^cc ^ 1/t In t, 
from an inefficient, memory-less search of all oscillators, 
for an ever-decreasing un-relaxed population. 

Figure 10 shows for this model, the very different ac¬ 
ceptance fractions Aacc{t) versus time t for the three 
Regions. At T = 0.76 in Region 4, Aacc{t) is nearly 
zero during incubation. The acceptance spikes at con¬ 
version times t = tjn during conversion from domain- 
wall vapour to liquid, and falls again to zero beyond 


t = tc, in the crystal phase. The spike occurs at the 
same time as the sharp rise of the martensite fraction 
through nm{tm) = 0.5, giving a physical justification to 
this earlier definition of 
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FIG. 10. Monte Carlo Aceeptanees : Plot of Aacdf) versus t 
for the three temper at ure-quench Regions 4, 2, 1. a) T = 0.76 
in Region 4; b) T = 0.55 in Region 2; c) T = 0.4 in Region 1. 


For T = 0.55 in Region 2, Aacc{t) is high initially, 
accepting most of the flips; decreases when the domain 
wall liquid phase is reached; peaks again near t = tc 
on domain wall orientations; and finally falls to zero ac¬ 
ceptances in the crystal phase. There can be a second 
peak just before t = tc^ where austenite droplets are 
generated to bind to the domain wallJ^. For T = 0.4 
in Region 1, the acceptance in the liquid phase is spiky, 
during domain wall motion in the liquid phase. Then 
there is a peak in Aacc{t)^ ^s a large number of austenite 
droplets or hotspots are generated to catalyse domain- 
wall symmetry-breaking orientational^, finally falling to 
zero as before, in the crystal. 

The ageing state has MC acceptance fractions that are 
nearly zero, and acceptance spiking at tm or tc is thus 
a diagnostic at all temperatures for the crossing of an 
entropy barrier. 


V. DOMAIN-WALL THERMODYNAMICS AND 

EFFECTIVE TEMPERATURES 

We approximate the MC Hamiltonian H (t) by that of 
independent spins at T, in a time-dependent local mean- 
fielcP^ cr(r, t), where the non uniformity comes from the 
domain walls. This implicitly assumes (consistent with 
simulations) that there is a separation of time scales, with 
individual spins flipping rapidly in response to quenches 
of the temperature, with domain wall configurations 
evolving more slowly. The Appendix obtains, within a 
‘time-dependent local mean-field’ approximation, expres¬ 
sions for the free energy E ELMF{t)^ internal energy 
I/(t), and entropy 5'entr(^) = —PLMF{t)-\-TS{t), in terms 
of the {5'(r, t)} configurations at a given MC sweep la¬ 
belled by t. 

We regard the domain-walls under stress as internal 
work sources^ that run freely after a quench. The overall 
change in the internal energy dU (t), by a First Law type 
relation, is a sum of contributions from the work done by 
domain walls (fLF(t), and the heat release dQ{t)\ 

dU{t) = d:w{t)pd:Q{t). ( 6 . 1 ) 












We need relations between the textural thermody¬ 
namics and increments in the heat and the work, at 
constant T. The heat release by spins, that are at 
the bath temperature, is taken as dQ{t) = TdSentr{t). 
For equilibrium, the Helmholtz free energy change be¬ 
tween thermodynamic states, is the available work at 
constant temperatur^^. We assume the free-running 
work increment at constant temperature saturates a sim¬ 
ilar availability, set by the evolving free energy change: 
dW{t) = dFLMF{t). At long times after entropy barri¬ 
ers are crossed, thermodynamic equilibration dF = 0 is 
accompanied by mechanical equilibration dW = 0. 

The evolving work rate is kF(t), and the heat emission 
rate is Q(t), where ‘rates’ are X{t 1) — X(t) = X{t). 
Fig 11 shows the probability distributions P{W), P{Q) 
over MC runs, for the rates of internal work done or heat 
emitted. They can peak at negative values, but finally 
both equilibrate to peaks centred at zero. Fig 12 shows 
that the mean rates of work done and heat emitted are 
suppressed in the ageing state by entropy barriers, but 
show large and sudden releases, as the entropy barriers 
are crossed. On equilibration, all mean rates tend to zero. 
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FIG. 11. Evolution of distributions for rate of heat emission 
and internal work : Plot of evolving P{Q,t) versus Q and 
P{W,t) versus W for quenches into the three Regions 4, 2, 
1. Top row: T = 0.76 in Region 4. Middle row: T = 0.55 in 
Region 2. Bottom row: T = 0.4 in Region 1. For all T, the 
final distributions are peaked symmetrically around zero. 



FIG. 12. Evolution of averaged rates of heat emission and 
internal work: Log-linear plots shows the mean rates <W > 
and < Q > versus the time t. For T = 0.76, the rates are 
zero in the ageing state, have large releases near tm{T), and 
finally vanish on equilibration. For T = 0.55,0.4 the curves 
are similar, but moved to the left, as conversion delays tm, 
vanish. For all T, the final mean rates are zero. 


tracked, by defining a Te//(t), through 
dW{t)^[l-X]dU{ty, dQ{t) ^ XdU{t). (6.2) 

-^eff J-eff 

With the previous increment relations this is manifestly 
equivalent to a ’microcanonicar definition, T/Teff(t) = 
TdSentr{t)/dU{t), as invoked in protein folding model 
The incremental work can also be related to changes 
in the local mean-field, 

(f ir(i) ~ - y] p(f, t)6a{r, t) (6.3) 

r 

where the local stress is p{f, t) = dFpMF{t)/d(T{r^ t). For 
static textures at equilibrium (7(r, t) = d(r) satisfying the 
mean-field self-consistency conditiorP^, the local stresses 
vanish in the domain-wall crystal phas^^, as seen in the 
Appendix. FiomdW ^ (1—T/Te//) of (6.2), this vanish¬ 
ing of trapped stresses is consistent with Teff(t) T. A 
detailed study of this final equilibration would involve the 
second entropy barrier at tc, and so the trapped-stress 
related effective temperature will be pursued elsewhere. 


VI. SUMMARY 


The energy of trapped stresses can only escape as heat 
released to the bath, since the boundary conditions are 
periodic, and not piston-like. One would like to relate 
a stress-induced heat release to re-equilibration of some 
effective temperature^ that has otherwise been defined in 
terms of the Fluctuation-Dissipation relatioiP^. 

In the equilibrium case, and for some externally im¬ 
posed work protocoP^, one can distinguislP^ between 
heat changes (that are occupancy changes of given en¬ 
ergy levels), and work changes (that are energy levels 
changes at fixed-occupancy). For temperature quenches, 
however, there is no such external sequential control, and 
both work and heat changes occur together. The relative 
proportion of spontaneous heat and work in ( 6 . 1 ) can be 


We develop a detailed understanding of the re¬ 
equilibration process of domain walls, in a martensite- 
related three-state model with powerlaw anisotropic in¬ 
teractions. There is a natural appearance of concepts 
borrowed from protein folding, of golf holes and funnels; 
and from oscillator relaxation models, of Monte Carlo ac¬ 
ceptance fractions. As found earlier, domain-wall phases 
after a quench, evolve from a domain-wall ’vapour’ to a 
’liquid’, and thence to a ‘crystal’ of oriented walls. There 
is a temperature regime where the martensite conver¬ 
sion delay t^ dominates the total delay. The evolution 
from a vapour-phase zero-energy droplet in zero-energy 
austenite to a negative-energy liquid of martensite do¬ 
main walls, is best understood in Fourier space. The 
droplet has an isotropic gaussian structure factor, peaked 
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at the Brillouin zone centre, over a butterfly-shaped, 
small golf hole with a negative energy funnel inside. The 
incubation delays come from a search for anisotropies at 
zero energy to roll into the golf hole. At low tempera¬ 
tures, the golf hole is large, and conversion from seeds 
occurs almost immediately. At temperatures above tran¬ 
sition, there is golf hole topology change, preventing the 
roll-in, and suppressing the conversion of austenite. 

In the ageing state, the MC acceptance fractions 
are negligible; both work and heat rates are zero; and 
trapped local stresses are held in place by entropy 
barriers. On crossing entropy barriers at the death 
of ageing, there are spikes in the MC acceptance, and 
sudden releases of the trapped stress and heat. On 
the resultant thermal and mechanical re-equilibration, 
acceptance fractions for spin-flips again vanish; work 
and heat rates are zero; and the oriented domain walls 
are stress-free, with a related effective temperature 
going to the bath temperature. The scenario may have 
relevance to other models of athermal re-equilibration 
after a quench, such as glassy and granular models. 

It is a pleasure to thank Shamik Gupta, Uwe Klem- 
radt, Stefano Ruffo, VSS Sastry and Sandro Scandolo for 
useful conversations. NS thanks the University Grants 
Commission, India for a Dr. D.S. Kothari postdoctoral 
fellowship. 


Appendix : Time-dependent local mean-field 
approximation : 

The uniform^ static mean-field approximation is famil¬ 
iar for ferromagnets, and antiferromagnets (where it is a 
staggered magnetisation). However a static, local mean- 
field can faithfully reproduce the domain-wall textures 
from simulations^. This can be generalizecP^ to an MC 
time-dependent local mean-field approach, to describe the 
evolving domain walls. 

The model Hamiltonian of (2.1) has the bath temper¬ 
ature T entering H(T) through the martensitic strain 
magnitudes e{T). With the partition function as Z = 
^{S} and the canonical free energy as F = 

—TlnZ, the entropy is Sentr = —dFjdT^ the internal 
energy defined by U = F TSentry is then not just the 
averaged hamiltonian, but is 

U =<H>- < TdH/dT > . (Al) 

With strain-pseudospin patterns {5'(r, t)} evolving un¬ 
der Monte Carlo dynamics, the weight factor e“^/^is 
truncated withirP^ a ‘time-dependent local mean-field 
approximation’ that we restate here for completeness. It 
is defined by the substitution into the coordinate space 
Hamiltonian as 

S{r)S{F) S{f)a{F,t)-\-(T{r,t)S{r)—(T{r,t)(T{r,t). {A2a) 


Here, the mean field spin a is defined as a run-average of 
a local spin variable at a site r, and an MC time t: 

(j{r,t) =< S{f,t) > • (A 26 ) 

The local mean-field weight is then e~^l^ 
where 

Hlmf/T = {A2c) 

r 

depends on individual spins in a mean-field, and 

= iA2d) 


where 

Qoir-f) = Do[{gL{T)+^oA/}Sf^^f^, + Fv{f-f^)]. {A2d) 

As mentioned in the text, the individual spins are as¬ 
sumed to respond instantaneously to the quenched tem¬ 
perature of the heat bath and to the influence of domain 
walls, that themselves evolve much more slowly, under 
MC dynamics. 

The corresponding substitution in the Fourier space 
Hamiltonian of (2.1) is 

\S{k)\‘^ 5'(^)cr(fc, t)*+cr(fc, t)5'(^, t)* —|cr(^, t)p. (A3a) 

Here 

Hlmf/T = J2(lik,trS{k,t)-lqo{k,t)W{k,t)f, {A3b) 

k 

where 

q{k, t) = qQ{k)a{k, t), {A3c) 

and 

qo{k) ^ Do[{gL{T) + eoK^} + Fv{k){l - {A3d) 

with Do{T) = Kq/T = 2e^{r)Eo/T. 

The thermodynamic functions are all taken as zero in 
F = 0 uniform austenite, and the approximate LMF free 
energy is 

Flmf = -T{y;ini[l+2coshg(r^]-i y;go|cr(fc)p}, (^4) 

r k 


The internal stress p{rA) = QElme/ d(7{rA) is 

p«t)=- E ,) 01 <-<5) 


and vanishes at the self-consistent, static equilibrium 
texture^^(j(r, t) = d(r), when the square bracket is zero. 

From (Al) and the Hamiltonian (2.1), the internal en¬ 
ergy is 
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u = [1- <H> E < >’ (^6) 

k 

where the averages are with the LMF weight. The en¬ 
tropy is taken as the difference of (A4), (A6) 

TSentr = —FlMF + U. {A7) 

With Td/dT = [r(T) — r(0)](i/(ir and the definitions of 
the text, Tde^/dT = — [r(T) — r(0)]/(4yT^^^37^); and 
TdgL/dT = [r(T) - t(0)][1 - 2{£2 _ l)/(4v/l - 3r/4)]. 

Of course, these expressions are in terms of the time- 
dependent local mean-field cr(r, t). We sidestep the evalu¬ 
ation through (A2b) of cr(r, t) at each time t, by invoking 


the spirit of mean-field approximations, namely that 
‘the function of an average is an average of the function’, 
and taking 

^)}) —< FLMF{{S{r,t)}) >; (ASa) 

l/({a(f, t)}) c^< U{{S{r, t)}) >; {A8b) 

where the averages are now taken over each distinct MC 
runs. 

Thus the (time-dependent) LMF expressions of (A4), 
(A6), (A7) yield expressions for the domain-wall ther¬ 
modynamics of an evolving texture {5'(r, t)} over a each 
distinct re-equilibration run, that can then be averaged 
over many runs. This approach has been used for Figs 
11 , 12 . 
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